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"Method for determining reflections in an area". 
* * * 

DESCRIPTION 

The present invention refers to a method for determining reflections in 
5 complex environments, by means of ray tracing 

The rapid tracing of rays, that follow reflection laws similar to 
geometric optics laws, is a problem of considerable interest in numerous 
fields of application. In particular, the techniques capable of carrying out 
this operation find immediate application in the simulation of acoustic 
10 reverberations in complex environments, in the photo-realistic rendering of 

models of three-dimensional scenes, or in the prediction of the 
electromagnetic propagation of wireless systems in the presence of multiple 
reflections. 

The acoustic or electromagnetic path curves follow paths that, with 
1 5 excellent approximation, meet the geometric optics laws. According to these 

laws, a ray that has an incidence on a planar reflector with angle a in relation 
to the normal to the plane is reflected by this in coplanar direction with the 
plane formed by the incidence vector and the vector normal to the reflector, 
which results specular in relation to the latter vector (output angle equal to a 
20 -a). A ray can undergo many reflections by the walls present in the 

environment, each time diminishing until it becomes of negligible entity. 
Given a source and a receiver however located in an environment with 
various reflectors, many paths can exist (in direct or indirect visibility) that 
connect these points. 

25 There are various approaches to the problem of tracing rays, such as the 

image sources method, the ray tracing method, and the beam tracing method. 

The image sources method calculates the paths corresponding to 
specular reflections considering the virtual sources obtained through the 
specular reflection of the sources in relation to all the polygons that 

30 represent the reflectors of the environment. For each virtual source a 
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specular reflection path can be constructed through the iterative intersection 
of a rectilinear segment that links source with receiver through the existing 
reflection planes. The specular reflection paths can be calculated up to an 
arbitrary order through the recursive generation of virtual sources. The main 

5 advantage of the image method is its solidity, as it guarantees that all the 

required paths are determined up to a pre-assigned reflection order. Against 
this, the cost of calculating the method grows in exponential proportion of 
the reflection order and the control of the visibility conditions (to be carried 
out for each of the virtual sources) results to be of limited complexity only 

10 for environments that are not too complex (in fact conditions of occlusion, 

plausible reflection, etc. have to be considered). 

The ray tracing methods are based on the direct tracing of the ray 
starting from the source until it reaches the receiver o vice versa. Whatever 
point the ray tracing starts from, these methods typically require the 

1 5 exploration of a finite, but high, number of output directions, to verify that a 

path exists that complies with the geometric optics laws and that reaches the 
arrival point after a certain number of reflections. However, this approach 
often results to be burdensome to calculate as in order to be certain not to 
excessively undersample the space of the possible directions to explore (and 

20 erroneously miss the arrival point), the number of tracing attempts has to be 

very high. 

Beam tracing is capable of preventing this problem through the tracing 
of compact beams of rays that originated from the same source (real or 
virtual). This eliminates the risk of missing the arrival point in the tracing 
25 operation. 

When a beam has incidence on a reflector, all the rays that compose it 
are reflected according to the geometric optics laws forming a reflected 
beam that results in coinciding with the specular version (in relation to the 
surface of the reflector) of the beam that would be obtained if the reflector 
30 was not present 
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In general, however, the beam has incidence only in part on the 
reflector, with the result that only a portion thereof is reflected while the 
remaining part proceeds undisturbed in its original path. In this case the 
beam separates into two or more sub-rays in correspondence with the 
5 reflector. To carry out beam tracing thus means to construct the tree of the 

reflections and of their branches that the beam undergoes in its path in space. 
Clearly, the ray undergoes a natural attenuation that grows with the increase 
of the length of the path followed. A further attenuation is introduced by the 
reflection walls. This enables the growth of the tree to be interrupted where 
10 the attenuation of the beams reaches a preset threshold. 

The reflection tree or beam tree is constructed so as to gather in 
gerarchic mode all the information needed for the reconstruction of the 
acoustic paths existing between the source and an arbitrary point. In 
particular, for each branch point of the tree the position of the source (real or 
15 virtual) that generates the reflection to that branch point has to be specified, 

the geometrical parameters that characterise the extension of the beam, and 
the information need to identity the reflector - if any - on which the beam 
will have incidence in following its path. It is assumed, without loss of 
generality, that the beam has been preliminarily divided into adjoining sub- 
20 beams. 

Once the beam tree associated to a given position of the source has been 
constructed, the determination of all the beams that hit a certain point of the 
space (receiver) is immediate. Going back up the tree it is then possible to 
derive the direct or reflected paths that link it to the source. This operation 
25 results, in fact, in being very rapid as it only requires consultation of the data 

structure used to memorise the beam tree. 

It can be seen however that a movement of the source or even of only 
one of the reflectors entails the recalculation of the tree, an operation that is 
generally more burdensome than its simple consultation. 
30 In view of the state of the technique described, object of the present 
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invention is to speed up the construction of the beam tree when the source 
moves, through the separation of the phases of analysis of the visibility and 
tracing of the beams. 

In its first aspect the present invention refers to a method for 
5 determining the region of visibility between at least a first reflector and a 

second reflector comprising the following phases: representing said first and 
second reflector in a system of coordinates (x, y, z); said method is 
characterised in that it further comprises the phases of: carrying out a similar 
transformation of said system of coordinates (x, y, z); detemiining the region 
10 of visibility of said second reflector in relation to said first reflector as the 

set of the parameters of the straight lines that link a generic point of said first 
reflector with a generic point of said second reflector. 

Preferably, the affine transformation comprises the execution of a 
rotatory translation of said system of coordinates (x, y, z) so that said first 
15 reflector is placed on the plane x = 0. 

Preferably, the affine transformation comprises the execution of a 
scaling down of said system of coordinates (x, y, z) so that said first reflector 
assumes a preset dimension. 

Advantageously, the method comprises a third reflector; and the 
20 following phases: determining the region of visibility of said second and 

third reflector seen from said first reflector; in the case that there are 
overlapping regions of visibility: tracing a semistraight line in said system 
of coordinates that has parameters comprised in the overlapping region 
starting from said first reflector; determining what is the succession of the 
25 intersections between said second and third reflector; assigning the 

overlapping portion of the region of visibility to the reflector having priority 
succession order; said regions of visibility represent a visibility diagram. 

In particular the method is characterised in that it represents said first 
reflector and second reflector in a system of coordinates (x, y) by means of 
30 segments; carrying out an affine transformation of said system of 
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coordinates (x, y) capable of bringing said first reflector to assuming 
coordinates of the extremes in the points (0, m) and (0, n) with m < n; 
representing said second reflector by means of the following system of 
equations x = et + f,y=gt + h, with 0 < t < 1 ; representing a generic 

5 straight line by means of the parameters a, b of the equation y = a x + b; 

determining the region of visibility by determining all the straight lines that 
pass through a generic point of the first and of the second reflector 
combining the previous equations and obtaining the following system g t + h 
= a (e t + f) + b, with 0 < t < 1 and with m < b < n. 

10 In its second aspect the present invention refers to a method for 

determining a beam tree of the beams of the rays on a plurality of reflectors 
that comprises the phases of prior determination of the set of the regions of 
visibility from each reflector in accordance with the previously described 
method. 

1 5 Advantageously, the method for determining a beam tree of the beams 

of the rays on a plurality of reflectors comprises the phases of prior 
determination of the visibility diagram in accordance with the previously 
described method; positioning a source in a system of coordinates (x, y); 
determining the areas of the reflectors illuminated by said source in said 

20 system of coordinates (x, y); memorising the coordinates of said areas of the 

reflectors illuminated by said source; representing a virtual source for each 
reflector Uluminated; applying said aflRne transformation of said system of 
coordinates (x, y) to said virtual source and to said illuminated region; said 
transformed illuminated region is represented by means of a segment of 

25 extremes (0, b 0 ) and (0, bi); representing said transformed virtual source in 

the space of the parameters (a, b) by means of the equation y = a x + b, and 
the illuminated region by means of the disequation bo < b < bi ; the system of 
equations y = ax + b and bo < b < bi represents a segment of illumination in 
the space of the parameters (a, b); intersecting said segment of illumination 

30 with said visibility diagram obtaining sub-segments and thus sub-intervals of 
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the interval bo < b < bi; said sub-intervals will represent portions of the 
illuminated region and each one will illuminate a new reflector. 

In its third aspect the present invention refers to a program for computer 
comprising a program code that carries out all the phases of the method for 
5 determining the region of visibility between at least a first reflector and/or of 

the method for determining a beam tree when the program is carried out on a 
computer. 

In its fourth aspect the present invention refers to a program for 
computer recorded on a support that can be used by a computer for 
10 controlling the execution of all the phases of the method for determining the 

region of visibility between at least a first reflector and/or of the method for 
determining a beam tree. 

The present invention permits to carry out the bulkiest calculation 
operations of assessment of the visibility in a preliminary phase that 
1 5 considers only the reciprocal positions of the reflectors, independently from 

the position of 1he source. The visibility is in fact redefined in a suitable 
parametric space, where the information for its assessment can be organised 
under the form of diagrams. The construction of the beam tree is then carried 
out each time the source moves, using the precalculated information of 
20 visibility. 

The characteristics and advantages of the present invention will appear 
evident from the following detailed description of an embodiment thereof, 
illustrated as non-limiting example in the enclosed drawings, in which: 
Figures la, lb and 1c show an example of rotatory-translation, of 
25 normalization and transformation from a space with coordinates (x, y) to a 

space with parameters (a, b); 

Figures 2a and 2b show an example of construction of the visibility 
diagram; 

Figures 3a and 3b show an example of construction of a level of the 
30 beam tree. 
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As the method herein proposed refers to the tracing of paths that follow 
the geometric optics laws, phenomenon such as diffraction will not be 
considered, which occurs when the dimensions of the reflectors become 
comparable with the wavelength of the wave that propagates in the 

5 environment. This assumption, even though correct in the case of high 

frequency electromagnetic waves (light or microwaves), results only 
approximately satisfied in the case of the propagation of acoustic waves, 
and in feet it has been calculated that these models are valid for simulations 
above 150Hz (corresponding to the wavelength of approximately 2 metres, a 

10 length that starts to be comparable with the dimensions of the walls of a 

living or working environment). In addition, even though it is true that 
sufficiently wide walls reflect in an almost specular manner, rough edges of 
rigid reflectors tend to generate a diffraction that is not negligible. In 
particular, small objects (a vase, a cup, etc.) tend to have significant effects 

15 on the acoustic field for sounds whose frequency is above 4kHz. These 

objects, however, are not generally considered in the modellation of acoustic 
environments. 

For simplicity of description, hereinafter the phenomenon of the 
transmission through the walls will also be neglected. Hereinafter it will be 
20 illustrated how it is possible to take account of it through simple 

modifications of the method illustrated herein. 

The entities that we shall now consider are the sources, the receiver and 
the reflectors. 

It is assumed that the sources and receivers are pointed elements, while 
25 the reflection walls are described by segments in the two-dimensional case 

and by portions of polygonal planes, for example triangles, in the three- 
dimensional case. Each wall of reflection defines two reflectors, one for 
each side. This distinction, as will be seen, is substantial in that it has 
consequences on the choice of parameters that characterise the visibility. 
30 We can say that a source "illuminates" a point if there is a direct path 
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that links them. With the term 'Virtual source" we intend the specular 
reflection of a real source in relation to a reflecting wall. In addition that 
region of a reflector that is directly illuminated by the source considered is 
defined "active". The active portion of a reflector can be composed of one or 
5 more active segments (also disconnected because of possible partial 

occlusions). Each of the active segments generates a new reflected beam. 

Identifying the reflectors present in the environment with a whole index 
i, the visibility is defined as a function that associates the index of the visible 
reflector to a view point and to a view direction. This function is thus 
10 defined in the space of the parameters that characterise view point and 

direction, and is constant at times. The visibility from a reflector is thus 
defined as the visibility function from all points of the reflector itself. 

The assessment of the visibility of an environment, that is of a set of 
reflectors, consists in the construction of the visibility functions associated 
15 to all the reflectors of the environment. As it is assumed that the reflecting 

walls reflect on both sides, the visibility of an environment composed of M 
walls will be characterised by 2M visibility functions. 

The preliminary assessment procedure of the visibility of the 
environment through the construction of the visibility by the generic 
20 reflector is now illustrated To be able to conveniently define the space of 

the parameters on which the visibility function assumes values, it is 
preferable to carry out a affine transformation of the geometric space that 
leads the reflector to assume preselected position, orientation and, finally, 
shape. 

25 As previously illustrated, the visibility function assumes values on the 

space of the parameters that characterise a view point and direction. If the 
point lies on the reflector, the space of the parameters will be four- 
dimensional in the case of three-dimensional environments, or two- 
dimensional in the case of two-dimensional environments. We assume that 

30 the geometric space has been preliminarily submitted to affine 
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transformation that brings the reflector to lie on plane x = 0, with reflecting 
face turned towards the semi-axis x > 0. A convenient parametizing for view 
point and direction is thus given by the system of equations 
y=ax+b 
5 z=cx+d 

where the couple of parameters (b, d) identifies the view point on the 
plane (y, z) of the reflector, while the couple of parameters (a, c) defines the 
view direction. In the case of two-dimensional environments only the first of 
the two equations has to be considered, and the parameters a and b will 
1 0 identify the view direction and point respectively. 

For reasons of simplicity of description, we will limit ourselves for the 
moment to consider only the case of two-dimensional environments. In this 
case, the affine transformation of the geometric space can be defined simply 
as a rotatory-translation with scaling down that transforms the reflector in 
1 5 the segment that links the points(0, -1) and (0, 1), with reflection face 

turned towards the semi-axis x >0. This enables the space of the parameters 
at the interval -1 < b < 1 to be limited The visibility function will thus be a 
constant function at times in the space (a, b), where the values identify the 
visible reflectors. 

20 We will consider two reflectors ro and r i, having as references 

respectively q 0 , po and qi, pi, located casually in a space with coordinates (x, 
y), as shown in Figure la. 

We roto-translate the reflectors placing the reflector ro in question to 
coincide with the vertical axis y (x=0). We normalize, again for practicality, 

25 the reflectors, that is the segment having as references p'o and q f 0 of 

coordinates (0, 1) and (0, -1) has been chosen as representation of the 
reflector r 0 . Thus a representation of the reflectors as shown in Figure lb is 
obtained. Representing the reflectors in the space of the parameters (a, b) the 
diagram in Figure lc is obtained 

30 In the space of the parameters (a, b), in Figure lc, the space o comprised 
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between the four straight lines p"o, p f, i, q"o, q"i represents the visibility of r ! ! 
from r'o- 

We define the region of visibility of a generic reflector in relation to the 
reflector of reference as that region of the space of the parameters (a, b) to 
5 which view points and directions from which that reflector is visible 

correspond. As occlusions by other reflectors can occur, this region can 
result empty, or composed of one or more convex polygons. In each region 
of visibility the visibility function will obviously result to be constant Vice 
versa, the region of visibility associated to the umpteenth reflector, 
10 corresponds to the region in which the function assumes value i. It is in 

addition useful to define the region of potential visibility of a reflector in 
comparison to the reference reflector, like the region of visibility calculated 
in absence of all the other reflectors. The effective region of visibility will be 
justly contained in the region of potential visibility, as the regions of 
15 potential visibility of all the reflectors can overlap. 

The construction of the visibility function associated to a reference 
reflector is carried out in two successive phases: first the regions of potential 
visibility of the single reflectors are assessed, then the regions of effective 
visibility are assessed through determinations of the reciprocal occlusions. 
20 For determining the region of potential visibility of the generic reflector 

compared to the reference reflector, the parameters of all the straight lines 
that start from any point of the reference reflector (in positive direction) and 
have incidence on the reflecting side of the reflector considered at any point 
are determined. To this object, first of all the generic reflector is defined 
25 with the couple of equations: 

x=et+f 
y-gt+h 

where t is a parameter comprised between 0 and 1 . Replacing these 
expressions in the equation y = ax + b we obtain a region in the space of the 
30 parameters characterised by: 
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gt + h = a(et + f)+b, 0<t<l 
-l<b<l 

This region results to be thus the intersection between a beam of semi- 
straight lines and the reference strip -1 < b < 1 . 

5 When the areas of potential visibility of two reflectors overlap each 

other in the space of the parameters, it means that a reflector partially blocks 
the other in relation to the generic view point on the reference reflector. The 
second phase of the procedure for assessing the visibility function consists in 
assigning the correct index to the regions in which the regions of potential 

10 visibility overlap each other. This operation can be carried out in various 

ways. The most simple method consists of determining which reflector 
occludes all the others by means of tracing a sample view ray chosen inside 
the overlapping region. At the end of this operation we obtain a partitioning 
of the space of the parameters in areas of effective visibility, to each of 

15 which is associated at the most a visible reflector. We will call this 

geometric representation of the visibility function, visibility diagram. 

Now we consider four reflectors ri-r4, and we want to determine the 
visibility diagram of the reflector r 4 . Thus, we roto-translate and normalize 
the reflector r 4 and we obtain a representation in space with coordinate (x, 

20 y) as shown in Figure 2a, and in the space of the parameters (a, b) as shown 

in Figure 2b. We also consider a source S represented in the space of the 
parameters (a, b) with the straight line S\ It can be seen that in Figure la, the 
reflector ti partially occludes the reflector n. An occlusion in the space of 
the parameters (a, b) is represented by an overlapping of areas. Analysing 

25 tiie slopes of the straight lines that enclose the overlapped areas it is possible 

to determine in the correct way the relative visibility. In particular, the 
greater slope dominates the minor slope. 

This operation of construction of the functions (diagrams) of visibility is 
carried out on all the reflectors present in the scene (two per wall) in a 

30 preliminary analysis phase. 
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At this point the real tracing of the beams for the construction of the 
beam tree can be started, an operation that is possible when the position of 
the source is known. As already partially anticipated, each branch point of 
the beam tree is associated to a beam, and is characterised by the following 
5 information: position of the real or virtual source, index and illuminated 

portion of the reflector on which the beam has incidence. 

The present method enables in particular the construction of the upper 
branching of the beam tree at the first level to be sped up. The construction 
of the branching of the root of the tree (primary beams) can come about, as 
10 in the classic method, tracing a beam at 360 degrees from the source and 

verifying the occlusions in the geometric space. 

The beam tracing procedure comes about in recursive mode. After 
having divided the beam emitted by the real source into primary beams, each 
associated to a reflector in visibility, they are reflected and the recursive 
15 determination of the successive branching is made. 

Given a beam reflected by the umpteenth reflector, the information of 
the visibility diagram associated to it is used to determine the subdivision of 
the beam into sub-beams, one for each reflector in visibility. 

In practice, the procedure is as follows: first of all the transformation 
20 used to bring the reflector associated to the beam in canonic reference 

position is applied to the virtual source, that brings it to lie on the left of the 
reflector. At this point the interval b 0 < b < bi that characterises the portion 
of beam that actually illuminates the reflector is determined. Then the 
representation of the beam on the dominium of the visibility diagram is 
25 identified. 

Replacing the point identified by the virtual source in the equation 
y = ax + b the equation of a straight line in the space of the parameters (a, b) 
is obtained, which together with the condition bo < b < bi represents a 
segment. Intersecting this segment with the areas that compose the visibility 
30 diagram of the reflector a series of sub-segments of the original segment is 
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obtained, each of which represents in the geometric space a new beam that 
will have incidence on a single reflector. Using this information a new level 
of the tree will be created. 

Figure 3 refers to the construction of a level of the tree, based on the 
5 reflections of the reflector r 4 of Figure 2. Positioning the straight line S 1 in 

the space of the parameters (a, b), this intersects the reference area (-1 < b < 
1) in 5 sub-rays numbered from 1 to 5. The rays 1 and 3 are rays that are not 
reflected and proceed without limits, while the rays 2, 3 and 5 are reflected 
by the reflectors r2, ri and rs respectively. The information relating to the 
10 position of the source that generates the reflection to that branch node, the 

geometric parameters that characterise the extension of the beam, and the 
information needed to identify the eventual reflector on which the beam will 
have incidence in continuing its path are associated to each branch node. In 
this manner the tree is created, and the process is repeated for each reflection 
1 5 of the rays, until the number of reflections reaches a preset number or the 

energy associate to the ray becomes negligible. 

The information present in the beam tree enable the generation of a 
filter (FIR) on the basis of which paths reach the location of the listener. 
This updating is rapid and can be carried out frequently. 
20 Through this process it is possible to obtain for each path its own 

direction of arrival (DOA), the length of the path and the list of the reflectors 
with which it has interacted along the path. This phase requires all the 
branches and the leaves of the tree to have been visited and has to be 
controlled if the receiver is present in the beam being considered. If the 
25 receiver is present, the direction of arrival (DOA) of the ray and the length 

of the path are calculated. Therefore, passing along the whole tree to the root 
the list of reflectors met is obtained. 

Since it is impossible to produce a flow for each DOA, angular intervals 
have been grouped together to generate a limited number of audio flow, one 
30 per interval. Each of these angular intervals will be attributed to separate 
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delay lines of the filter FIR, that are calculated taking into account all the 
paths and the relative attenuations. 

The filter calculated in this manner is capable, given an anechoic 
source (that is a source recorded in a room without reverberation), of 

5 obtaining the result of listening to that sound in a room represented by the 

geometry of the reflectors. 

The method allows immediate extension to the three-dimensional case. 
Also in the 3D case it is possible to cany out an affine transformation of the 
geometric space that enables the scansion interval to be preset, and to 

10 simplify the representation of the visibility function. As we will see 

hereinafter, without losing generalities, we will consider a geometric 
transformation that brings the plane on which the reflector lies to coincide 
with the plane x=0 and brings the polygon to assume a preset canonic form. 
As already previously mentioned, the space of the parameters of the 

15 straight lines in three-dimensional environments is four-dimensional. 

Downstream from the geometric transformation of the polygonal reference 
reflector, we can express the generic straight line of the space 3D with the 
system of equations y = a x + b and z = c x -4- d, in which (a, b, c, d) identify 
the space of the parameters. In this space the representation of the potential 

20 visibility is still given by set of the straight lines that have incidence on both 

the reflectors considered, which this time is a closed portion of the four 
dimensions space. For practicality, it is also possible to limit before the 
category of polygons that interest A particularly advantageous choice is that 
of the triangles: this transformation does not remove generalities from the 

25 method, as it is always possible to represent a complex plane polygon by 

means of its triangulation; in addition it is quite frequent, also in Computer 
Graphics, that the surfaces in the three-dimensional space are modelled with 
triangle nets. In the case of two triangular reflectors, the region of potential 
visibility can be expressed in the form of 6 disequation system, needed for 

30 the set of the straight lines having incidence on both the triangles. For 
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example, expressing the generic reflector by means of the system of 
equations: 

x = a x t + bxS 4c x 

y = a y t + byS +Cy 
5 z = azt+b z s+c z 

t>0 

s>0 

t + s< 1 

The system that represents the generic straight line passing through the 
10 reference reflector, previously normalized by an affine transformation that 

makes its three vertexes coincide with the three points (0, 0, 0), (0, 1, 0) and 
(0, 0, 1), is represented by the system: 
y=ax+b 
z=cx+d 
15 b>0 
d>0 
b + d<l 

Replacing the first three equations of the first system, in the equations 
of the second system we obtain two equations, from which we can derive the 

20 parameters s and t in function of the variables a, b, c, d. Replacing the values 

of the parameters s and t in the disequations of the first system, and adding 
the disequations of the second system we obtain a system of 6 disequations 
in the parameters a, b, c, d, that delimit the area of potential visibility 
between the two reflectors. 

25 In a similar manner to the two-dimensional case it is possible to decide 

the order of occlusion between overlapped regions of potential visibility, and 
thus it is possible to construct for each reflector a visibility diagram 
represented by a subdivision of the space of the parameters in regions of 
effective visibility, each of which is associated at the most to a reflector. 

30 The representation of the source and of the illuminated area is given by 
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describing all the straight lines that intersect it with a limited beam, 
obtaining a system of two equations and some dis equations in the four 
unknown quantities, that identifies a portion of a two-dimensional subspace. 
The subdivision of the beam will come about in this case intersecting the 

5 areas of the visibility diagram associated to the reflector that has generated 

the beam with the portion of plane that represents the beam reflected in the 
space (a, b, c, d). The method for the construction of the tree remains thus 
completely similar to the two-dimensional case. 

It can be noted that the choice of the transformation on the plane x = 0, 

10 both in the two-dimensional case and in the three-dimensional case, is not 

binding, but it has been introduced with the sole object of simplifying the 
exposition, placing the equations of the straight lines in an easier form. 

The method for determining the visibility between reflectors and the 
method for determining a beam tree, herein described, can be transcribed in 

15 a program, in a well-known manner by a technician of the sector, and can be 

memorised on any type of memory or support (floppy, CD) and/or can be 
made to be carried out by a computer. 



